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A note on fluctuations 
for internal diffusion limited aggregation" 

Amine Asselah 1 Alexandre Gaudilliere 2 

Abstract 



We consider a cluster growth model on Z d , called internal diffusion limited aggre- 
£NJ . gation (internal DLA). In this model, random walks start at the origin, one at a time, 

and stop moving when reaching a site not occupied by previous walks. It is known that 
p/ j the asymptotic shape of the cluster is spherical. Also, when dimension is 2 or more, 

and when the cluster has volume n d , it is known that fluctuations of the radius are at 
most of order n 1 / 3 . We improve this estimate to dimension 3 or more. In 

so doing, we introduce a closely related cluster growth model, that we call the flashing 
process, whose fluctuations are controlled easily and accurately. This process is coupled 
to internal DLA to yield the desired bound. Part of our proof adapts the approach 
of Lawler, Bramson and Griffeath, on another space scale, and uses a sharp estimate 
(written by Blachere in our Appendix) on the expected time spent by a random walk 
, inside an annulus. 
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1 Introduction 



x 

The internal DLA cluster of volume N, say A(N), is obtained inductively as follows. Initially, 
we assume that the explored region is empty, that is A(0) = 0. Then, consider N independent 
discrete-time random walks Si, ... , Sn starting from 0. Assume A(k — 1) is obtained, and 
define 

r fc = inf{t>0: S k (t) & A(h - 1)} , and A(k) = A(k - 1) U {S k (r k )}. (1.1) 

In such a particle system, we call explorers the particles. We say that the k-th explorer 
is settled on S^r^,) after time T k , and is unsettled before time tj.. The cluster A(N) is the 
positions of the N settled explorers. We study the growth of A(N), as N tends to infinity. 
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The mathematical model of internal DLA was introduced first in the chemical physics 
literature by Meakin and Deutch [llj. There are many industrial processes that look like 
internal DLA (see the nice review paper [5]). The most important seems to be electropol- 
ishing, defined as the improvement of surface finish of a metal effected my making it anodic 
in an appropriate solution. There are actually two distinct industrial processes (i) anodic 
levelling or smoothing which corresponds to the elimination of surface roughness of height 
larger than 1 micron, and (ii) anodic brightening which refers to elimination of surface defects 
which are protruding by less than 1 micron. The latter phenomenon requires an understand- 
ing of atom removal from a crystal lattice. It was noted in [11] , at a qualitative level, that 
the model produces smooth clusters, and the authors wrote "it is also of some fundamental 
significance to know just how smooth a surface formed by diffusion limited processes may 
be". 

Diaconis and Fulton [2] introduced internal DLA in mathematics. Their model is more 
general than ours: explorers can start on distinct sites, and the explored region at time 
is not necessarily empty. They were interested in defining a random growth process by 
iterating simple operation. They introduced many variations, and treat, among other things, 
the special one dimensional case. 

In dimensions two and more, Lawler, Bramson and Griffeath [5] prove that in order to 
cover, without holes, a sphere of radius n, we need about the number of sites of Z d contained 
in this sphere. In other words, the asymptotic shape of the cluster is a sphere. Then, Lawler 
in [7] shows subdiffusive fluctuations. The latter result is formulated in terms of inner and 
outer errors, which we now introduce with some notation. We denote with || • || the euclidean 
norm on M. d . For any x in R d and r in 1, set 

B(x,r) = {y G R d : \\y - x\\ < r) and M(x,r) = B(x,r) f]Z d . (1.2) 

For A C Z d , |A| denotes the number of sites in A. The inner error Si(n) is such that 

n - S z (n) = sup {r > : 1(0, r) c A(|1(0, n)\)} . (1.3) 

Also, the outer error 6o{n) is such that 

n + S (n) = inf {r > : A(|1(0, n)\) C 1(0, r)} . (1.4) 

The main result of [7] reads as follows. 

Theorem 1.1 [Lawler] Assume d > 2. With probability 1, 

5 An) $o{ n ) 

lim ~uj\ — r~Y2 = °> and hm ~rj3i — r\4 = °- L5 

n->oo n L ' 6 log(n) n->oo n 1 '^ log(n) 4 

Since Lawler's paper, published 15 years ago, no improvement of these estimates was achi- 
eved, but it is believed that fluctuations are on a much smaller scale than n 1 ^ 3 . Computer 
simulations [T2], [3] suggest indeed that fluctuations are logarithmic. In addition, Levine and 
Peres studied a deterministic analogue of internal DLA, the rotor-router model, introduced 
by J.Propp [1]. They bound, in [TU], the inner error 5j(n) by log(n), and the outer error 
5 (n) by n 1 " 1 ^. 

We present here an improvement on (11.51) in dimension d > 3. 
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Theorem 1.2 Assume d > 3. There is a positive constant Ad such that, with probability 1, 

8j(n) , , , Sn(n) 
lim sup v — < ^ lim _jj = 0. (1.6) 

n->oo log(n) n->oo n d+i log 2 (n) 

Let us now recall the approach of Lawler, Bramson and Griffeath in j8], explain our idea, 
and introduce a new growth model. The approach of [8] is based on estimating the number 
W(z) of explorers that visit each site z e Z d . It is based on the following observations, (i) 
If explorers would not settle, they would just be independent random walks; (ii) exactly one 
explorer occupies each site of the cluster. Thus, if we launch one explorer from each site of 
the cluster A(N), and call M(A(N),z) the number of crossings of site z, then M(NSq,z) 
which we define as W(z) + M(A(N), z) would be equal in law to the number of walks crossing 
z out of N independent walks started on 0. Even though M(A(N),z) and M(N5o,z) are 
dependent variables, some estimates on P(W(z) = 0) can be extracted from estimating the 
means of M(A(N), z) and of M(N5 , z). 

Rather than thinking in terms of one single site, we observe that a site has good chances 
to lie inside the cluster if some large region, about this site, is crossed by many explorers. 
How large should be this region, say C, and how many should be these crossings, say W(C), 
can be answered as follows. The size and location of C should be such that (i) the expected 
number of crossings of C is much larger than its standard deviation, and (ii) the number of 
crossings of C needed to cover C is of order \C\ (as suggested by the spherical shape result 

of m 

(i) E[W{C)\ > v / var(H/(C)), (ii) E[W(C)] > \C\. (1.7) 

Now, assume that |B(0, n)| explorers start at the origin. For a space-scale h(n) and an 
integer k > 1, to be determined, assume that B(0, n — kh(n)) is covered by settled explorers. 

Partition the shell S = 1(0, n - (Jfe - l)/i(n))\B(0, n - kh(n)) into about (n/h(n)) d - 1 
cells, each of volume h(n) d . Cells are brick-like domain, of side length the width of the 
shell. It is convenient to imagine that a cell C G S is the basis of a column of k cells 
reaching the boundary of B(0,n). It is also convenient to stop the explorers as they reach 
the boundary of B(0, n — kh(n)). Thus, with such a stopped process, explorers are either 
settled inside the B(0,n — kh{n)) or unsettled but stopped on its boundary, that we denote 
by <9B(0, n — kh(n)). What we have called earlier the number of explorers crossing C is taken 
here to be the unsettled explorers stopped on C fl dM(0,n — kh{n)). In these heuristics, we 
make the simplifying assumption that only explorers stopped on C fl dM(0,n — kh{n)) can 
cover C once released. 

• On the average, kh(n) d explorers are stopped on each cell of S. Thus, for k large 
enough (ii) of ( II. 7p should be fulfilled. 

• The standard deviation of W(C) is a delicate issue not yet settled. Here follows a heuris- 
tic justification for an upper bound for the standard deviation. If each explorer had to 
choose uniformly at random a cell C of S and were constrained to perform a reflected 
random walk inside a cone issued from the origin and with base C, then we would 
obtain another growth model whose fluctuations are expectedly larger than those of 
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internal DLA. For such a model the standard deviation of the crossing of a cell is of 



where r is the exit time from B(0, n — kh(n)), and Pq is the law of a simple random 
walk, S, starting at the origin. We expect f ll.Sp to be an upper bound for the standard 
deviation of W{C). 

At a heuristic level, (i) of (jl.7p follows if 



This discussion motivates a growth model associated from the start with the exponent 
■^j-g. We call this model the flashing internal DLA process, or simply the flashing process. 
The flashing process looks like internal DLA on a large scale, but has a distinct covering 
mechanism, which makes it much simpler to analyze. To obtain this growth model, we 
generalize the rules described in (11. ip . by enabling explorers to settle only at special times, 
called flashing times. Thus, each explorer is associated with a sequence of stopping times, 
and it is only at these times that it settles if outside the cluster. The precise definition of 
the chosen stopping times requires additional notation, which we postpone to Section El We 
describe here key features of the flashing process. 

First, 7L d is partitioned into concentric shells around the origin: a shell at a distance r 
from the origin has a width of order r^+i. Each shell is in turn partitioned into cells, which 
are brick-like domain, of side length the width of the shell. The key features are as follows. 

"Pa. An explorer flashes at most once in each shell. 

Vb. The flashing position in a shell, is essentially uniform over the cell an explorer first hits 
upon entering the shell. 

Vc. When an explorer leaves a shell, it cannot afterward flash in it. 

Feature Vb is the seed of a deep difference with internal DLA. The mechanism of covering a 
cell, for the flashing process, is very much the same as completing an album in the classical 
coupon- collector process. Thus, we need of the order of Vlog(V) explorers to cover a cell 
of volume V. For internal DLA, with explorers started at the origin, we need only of order 
V explorers to cover a sphere of volume V as shown in j8], and we believe that we need a 
number of explorers of order \C\ to cover a cell C. 

Feature Vc is essential for having the following coupling between flashing and internal 
DLA processes. 

Theorem 1.3 There is a coupling between the two processes such that, for all k, N > 1, and 
for a sequence {r^, fceN} going to infinity (that we describe in Section^), and h k = r^^ 1 ^ 



order 




(1.8) 




(1.9) 



• if A*(N) C 1(0, r fc + h k ), then A(N) C 1(0, r k + h k ) 
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• i/B(0,r fe + h k ) C A*(N), then B(0,r fc + h k ) C A(N). 



For the flashing process, we control easily the inner error. Then, to control the outer error 
we follow the approach of [7] , though with a simpler proof (allowed by the cell structure used 
to build our growth model). 

Theorem 1.4 There is a positive constant Ad such that, U.6]) hold for the flashing process. 
Also, we know that flashing internal DLA does exhibit power-law fluctuations. 
Theorem 1.5 With probability 1, 



Theorem 11.31 and Theorem 11.41 imply Theorem 11.21 

Let us now describe the heuristics behind Theorem 1 1.5 1 It is useful to organize the flow of 
explorers in the flashing process into exploration waves, in the way of Section 3 of [7] . That is, 
in the /c-th exploration wave, the explorers either stop as they reach the bulk of St, or settle 
before reaching Sk- Consider now the exploration wave associated with the last shell making 
B(0, n), say S*. Assume that at this time, the cluster fills B(0,n)\<S*. From our definition, 
S* has width of order n^, and the last shell receives a number of stopped explorers equal to 
its volume. There is necessarily one cell in S* which receives of the order of its own volume, 
and for a coupon-collector process, it is very unlikely that the explorers stopped in this very 
cell can cover the bulk of this cell before escaping S*. By feature Vc, if a hole is left in S* 
after the explorers leave S*, this hole remains uncovered forever. These heuristics pose two 
questions concerning internal DLA, we are unable to answer at the moment. 

• How many explorers stopped in the bulk of a cell are needed to cover the whole cell? 

• What is the correct order of fluctuations in internal DLA? 

The rest of the paper is organized as follows. Section [2] introduces the main notation, 
and recall well known useful facts. In Section [3j we build the flashing process, give an 
alternative construction, and prove Theorem 11.31 In Section HI we obtain a sharp estimate 
on the expected number of explorers crossing a given cell, and prove V-b. Both proofs are 
based on classical potential theory estimates. In Section [5], we prove Theorems 11.41 and 11.51 
Finally, in the Appendix, Sebastien Blachere gives a sharp estimate on the expected time 
spent in an annulus by a random walk. 

2 Notation and useful tools 
2.1 Notation 

We say that z,z' G Z d are nearest neighbors when \\z — z'\\ = 1, and we write z ~ z'. For 
any subset A C Z d , we define 



(1.10) 



0. 



dA = {z e Z d \A : 3z' eA,z'~z}. 



(2.1) 
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For any r < R we define the annulus 

A(r,R) = B(0,R)\B(0,r) and A(r, R) = A(r, R) n Z d (2.2) 

A trajectory 7 is a discrete nearest-neighbor path on Z d . That is 7 : N — > Z d with 7(2) ~ 
7(t + 1) for all i. The law of the simple random walk started in z, is denoted with P 2 . For 
a subset A in Z d , and a trajectory 7, we define the hitting time of A as 

H(A; 7) = min{t > : 7(f) G A}. 

We often omit 7 in the notation when no confusion is possible. We use the shorthand 
notation 

B n = B(0,n), B n = B(0,n), H R = H(B C R ), and H z = H({z}). 

For any a, b in R we write a Ab = min{a, 6}, and a V 6 = max{a,6}. Let T be a finite 
collection of trajectories on Z d . For _R > 0, z in Z d and A a subset of Z d , we call M(r, i?, z) 
(resp. M(r, i2, A)) the number of trajectories which exit B(0, R) on z (resp. in A): 

M(T,R,z) = J2Mah r )=z}, and M(r, i2, A) = J] M(I\ i2, z). (2.3) 

When we deal with a collection of independent random trajectories, we rather specify its 
initial configuration 77 e so that M(rj,R,z) is the number of random walks starting 
from 77 and hitting B(0, R) c on z. Two types of initial configurations are important here: (i) 
the configuration nl z * formed by n walkers starting on a given site z*, (ii) for A C Z d , the 
configuration 1a that we simply identify with A. For any configuration 77 G we write 

M = $>(*)■ ( 2 - 4 ) 



We are in dimension 3 or more, and Green's function of the simple random walk is well 
defined and denoted G. That is, for any x, y G Z d 



G(x,y) =E 3 



5^ l{S(n)=i/} 
n>0 



(2.5) 



For any A C Z d , we define Green's function restricted to A, Ga, as follows. For x, y G A 



G A (x,y) = E 2 



5^ l{5(n)=y} 
0<n<i?(A c ) 



(2.6) 



2.2 Some useful tools 

We recall here some well known facts. Some of them are proved for the reader's convenience. 
This section can be skipped at a first reading. 
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In [8] , the authors emphasized the fact that the spherical limiting shape of internal DLA 
was intimately linked to strong isotropy properties of Green's function. This isotropy is 
expressed by the following asymptotics (Theorem 4.3.1 of |9J). In d > 3, there is a constant 
Kg, such that for any z ^ 0, 

IGM-^I^ with C^^—, (2.7) 

where Vd stands for the volume of the euclidean unit ball in M. d . The first order expansion 
(12. 7p is proved in [9] for general symmetric walks with finite d + 3 moments and vanishing 
third moment. All the estimates we use are eventually based on (12.71) and we emphasize the 
fact that the estimate is uniform in ||z||. 

The following lemma is also used in the Appendix. 
Lemma 2.1 Each z* in 7h d \ {0} has a nearest-neighbor z (i.e. z* ~ z) such that 

\\z\\ < \\z*\\ -=. (2.8) 

2Vd 

Proof. Without loss of generality we can assume that all the coordinates of z* are non- 
negative. Let us denote by b the maximum of these coordinates and note that 

||-2*|| 2 < dfe 2 , and b > 1. (2.9) 

Denote by z the nearest-neighbor obtained from z* by decreasing by one unit a maximum 
coordinate. Using (12. 9p 

||^|| 2 - \\ z f = b 2 - (b- l) 2 = 2b- 1 > b > JO. (2.10) 

yd 

Note that <^M> follows from 2||z*||(||z*|| - \\z\\) > \\z*\\ 2 - \\z\\ 2 , and (12~T0|) . I 

We recall a rough but useful result about the exit site distribution from a sphere. This 
is Lemma 1.7.4 of [6]. 

Lemma 2.2 There are two positive constants c\,C2 such that for any z e dB(0,n), and 
n > 

Cl T <P (^ n )=^)<^ T . (2.11) 



n d x n 



Finally, we recall a well known large deviations estimate for independent Bernoulli vari- 
ables (see for instance Lemma 4.3 of [1]). 

Lemma 2.3 For any positive integer n, and {X%, . . . ,X n } n independent Bernoulli vari- 
ables, we have for any x > 0, and with X = X\ + • • • + X n 

max (P (X - E[X] >x),P(X- E[X] < -x)) < exp (- (mm (^j^j, |) ) ) ■ (2.12) 
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Remark 2.4 Note that for a sum of Bernoulli, var[X] < i?[X], and the following inequality 
is useful 



2 

X X 



max (P (X - E[X] >x),P(X- E[X] < -x)) < exp I - (^min (^^, 5 ) ) J ' (2 ' 13) 

Also, note that if E[X\ < E\Y\, where Y is a sum of m independent Bernoulli variables, 
then 



max (P (X - E[X] >x),P(Y- E[Y] < -x)) < 



3 The flashing process 



ex >(-( min (w)))' (2 - 14) 



In this section, we construct the flashing process. We then present a useful alternative 
construction of the same process. Finally, we prove Theorem 11.31 which couples the two 
processes. 

3.1 Construction of the process 

Partitioning the lattice. We partition the lattice into shells (<Sj : j > 0). For a given 
parameter ho > the first shell <5>o is the ball B(0, ho). The next shells are the annuli 

Sj = A(rj — hj,rj + hj), j > 1, (3.1) 

where rj and hj are defined inductively by r% — h\ — ho, and for j > 1 

1 

Tj + i — hj +1 = Tj + hj, and hj = r^ +1 . (3.2) 
We omit the easy check that (13. 2p yields 

We also define 

S = {0} and S 3 - = dM(0, rj ), j > 1. (3.4) 

Flashing times. Consider {Xj,Yj,j > 0} a sequence of independent Bernoulli variables 
such that 

P(Xj = 1) = 1 - P(Xj = 0) = 1, (3.5) 

P(Yj = 1) = 1 - P(Yj = 0) = j \ * j . = £ (3.6) 

Consider also a sequence of continuous independent variables {Rj,j > 0} each of which has 
density gj : [0, hj] — >• R + with 
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For j > 0, and Zj in Sj, let S be a random walk starting in Zj, an define a stopping time a 
as follows. If Rj = h for some h < hj then 



if Xj 



a 



{z,hA(rj + hj- \\zj\\)) c ) if X d = and Y d = 1, (3.8) 
H(A(rj - h, r j + /i) c ) if Xj = and Y} = 0. 



We set i7j = H(Ej), and we define the stopping times (<x,- : j > 0) as 

aj^Hj+aiSoe^), (3.9) 

where (9 stands for the usual time-shift operator. For j > we note that, by construction, 
S(t) G Sj for all t such that Hj < t < a and we say that <7j is a flashing time when S(<jj) 
is contained in the intersection between <Sj and the cone with base B(S(Hj), hj/2). We call 
such an intersection a cell centered at S(Hj), that we denote C(S(Hj)). In other words, for 
any z G 

C(z) = Sj f]{xeR d : 3A > 0, 3y e B(z, hj/2),x = Xy} . (3.10) 



The uniform hitting property The main property of the hitting time a constructed 
above is the following proposition, which yields property Vb of the flashing process to be 
defined soon. 

Proposition 3.1 There are two positive constants a.\ < a 2 , such that, for j > 0, Zj G 
and z* G C(zj). 

^<F Zj (S(a) = z*)<^. (3.11) 
The proof of Proposition 13.11 is given in Section |H 



The flashing process. Consider a family of N independent random walks (S* : 1 < i < 
N) with their hitting times, and stopping times {H it j, Zij, a^j : j > 0). Let also = S^H^j) 
be the first hitting position of on Sj. 

We define the cluster inductively. Set ^4*(0) = 0. For i > 1, we define r* as the first 
flashing time associated with S* when the explorer stands outside A*(i — 1). In other words, 

r* = min {<T itj : j > 0, S*(a id ) G C(z id ) nA*(i- l) c } , (3.12) 

and 

A*(i) = A*(t - 1) U {S*(t*)} . (3.13) 



3.2 Exploration Waves 

Rather than building A*(N) following the whole journey of one explorer after another, 
we can build A*(N) as an increasing union of clusters formed by stopping explorers on 
successive shells. Similar wave constructions are introduced in [8] and [7], with an equality 
in law between alternative constructions. However, the features of the flashing process are 
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such that in our case the two constructions are strictly equivalent. We use this alternative 
construction in the proof of Theorem 11.31 

We denote by G {X d ) N the explorers positions after the fe-th wave. We denote by 
A* k (N) and the set of sites where settled explorers are after the k-th wave. Our construction 
will be such that 

£ E fe & e u j<k Sj e> e A* k (N). (3.14) 

For k = we set f (i) = 0, and A* {i) = 0, for 1 < i < N. Then, for all k > 0, we set 
A* k+1 (0) = A* k (N). For i in {1, • • - , N}, we set the following. 

• If £ k (i) ^ Sfe, then 

6+i (*) = &(*) e U^S,-, and = *4fc+i(« - 1). 

• If f fc (i) G E fc and Si(a itk ) G C(z;, fe ) n ^(i - l) c , then 

fk+i(i) = Si(a i)k ) G «S fc , and = A* k+1 {i - 1) U {Si(a i>k )} . 

• If &(z) G E fc and $(0*,*) £ C(z i>k ) n ^(i - l) c , then 

6+1 (0 = Si(H i)k+ i) G Efc+i, and = *4£ +1 (z ~ !)■ 

In words, for each k > 1, during the fc-th wave of exploration, the unsettled explorers move 
one after the other in the order of their labels until either settling in S k -i, or reaching E& 
where they stop. We then define A*(N) by 

A*(N) = \jA* k (N). (3.15) 

fc>i 

We explain now why this construction yields the same cluster as our previous definition. 
An explorer cannot settle inside a shell it has left, and thus cannot settle in any shell Sj 
with j < k if it reaches Ti k . Now, since each wave of exploration is organized according to 
the label ordering, the fact that an explorer has to wait for the following explorers before 
proceeding its journey beyond E& does not interfere with the site where it eventually settles. 

3.3 Coupling internal DLA and flashing processes 

We use here the first definition of the flashing process, and realize the internal DLA process 
using the same randomness. 

Proposition 3.2 There is a coupling between the flashing and original internal DLA pro- 
cesses such that, for all N > 1, 

N N 

U {Si(t) : < t < Ti } C |J {S* (t) : < t < r*} (3.16) 

8=1 1=1 

and there is a one to one mapping ipx : A(N) — > A*(N) such that for z G A(N) 

if for k>l, z#\JSj, then ip N {z) £ |J Sj. (3.17) 

j<k j<k 
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Theorem 1 1.3 1 is a simple consequence of Proposition l3~2l On the one hand, if A*(N) C Uj^Sj 
for some k > 1, then, recalling that a flashing explorer cannot settle a shell it has left, the 
orbits of the N flashing explorers are all contained in Uj^Sj and, by (13.161) . so is A(N). On 
the other hand, if U,-<fc«Sj C A*(N), then, (13 . 1 7f) implying that with such a coupling 

\A{N) n U i<fc 5,-| > |.A*(iV) n U i<fc 5,-| = |U i<Jfe 5,-| , (3.18) 

which implies that Uj^Sj C A(N). 

Proof of Proposition [3721 We build the coupling together with the map ipN by induction 
on N. We use the trajectories of the flashing explorers to drive the internal DLA trajectories. 
We need a little more notation to do so. For each i < N, set for simplicity g*(i) = S*(r*), 
and denote by the length of the flashing trajectory (S*(t) : < t < t^) used to form 
the trajectories of the original explorers. Necessarily, we need < t* to have (I3.16p . For 
convenience, we partition A(N) into blue sites, say B(N), and red sites, say 1Z(N). 

We build a one to one map : A(N) — >■ {1, • • • , N}, together with the blue- red partition 
as follows. For each z in A(N), there are two possibilities. Either there is i < N such that 
Si( T i) = z an d ti,N = t*, and we say that z G B(N) and we set fisr(z) = i. Otherwise 
z G TZ(N). We then define Jn{z) as the label i of the flashing explorer that was driving the 
random walk Sj when the j-th explorer settled in z, and this will imply, by induction, that 

ti,N < r i- 

Finally the one to one map ipN '■ A(N) — >■ ^4*(A^) is the composition g* o f N . Note, first, 
that for all z G B(N), ip N (z) = z, second, B(N) C A(N) nA*(N) and last, 

z G B(N) & t fN[z)>N = r* fN(z) , z G K(N) & t fN{z)jN < T* fN[z) . 

We now start our induction with N = 1. The first trajectory {S{(t) : < t < t*) ends 
in g"*(l). We use this trajectory to build that of the first internal DLA-explorer. We set 
S'i(O) = 5*^(0) and immediately stop here since the origin = S'i(O) was initially unoccupied. 
As a consequence = 0, n = and there are two possibilities: either r* = 0, B(l) = {0}, 
K(l) = and /i(0) = 1, or, r* > 0, TZ{1) = {0}, 8(1) = and /i(0) = 1. 

Assume now that we have built from the trajectories {{S*(t) : < t < r*), % < A^}, 
the clusters *4*(iV), and the sets B(iV) and TZ(N), together with the times {i^jv, « < -/V}, 
and the one to one map /jy : ^4(iV) — >■ {1, • • • , iV}. We launch a new flashing explorer with 
trajectory (S^ +1 (t) : < t < t^ +1 ) that ends in g*(N + 1), and we start to define the 
(N + l)-th trajectory for the original internal DLA process by following S^ +1 : 

S N+1 (0) = S* N+1 (0), S N+1 (l) = S* N+1 {1), . . . (3.19) 

• If {£jv + i(£) : < t < r^r +1 } is not contained in -A(iV) then ^at+i settles the first time 
Sat+i exits A(N), that is at time (resp. on a site z) 

t N +i,N+i = mf{A; > : S* N+1 £ A(N)}, (resp. z = S* N+1 (t N+l>N+1 ) G A(N) C ) . 

We then set 

Vz < N, t itN+1 = t i>N , f N +iU(N) = fN, and f N +i{z) = N + 1, 
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and, if t N+hN+1 = r^ +1 (resp. t N+ i }N+1 < r^ +1 ). 

B(N + 1) ={z} U B(N) (resp. B(N)), 
K(N + 1) =K(N) (resp. {z} U TZ(N)), 

• If {S5r +1 (t) : < t < r^ +1 } is contained in A(N) = B(N) U K(N) then settles 
necessarily in a red site z since B(N) C *4*(iV). This red site is occupied by an explorer 
that was driven by a flashing explorer i = (z) when it settled, and we have < r* . 
After reaching Sjv+i^jv+i) = S'^ +1 (r^ , x ) in our definition of Sn+i, we set 

SW+i(l + ^+i) = S*{1 + t hN ), S N+1 {2 + r* N+1 ) = S*{2 + t hN ), . . . (3.20) 

we set t N+1>N+ i = rjy +1 , we turn blue the site z, we set f^ + i(z) = N + 1 and we 
proceed as previously: If {5*(t) : t^N < t < r*} <f_ A(N) then SV+i settles at the first 
exit from A(N) in some z' = S*(t it N + i) G A(N) C with U t N+i < t* and we set 

Vj < iV, j ^ i, tj jN+1 —tj jN , fN+i\A(N)\{z} — /ivU(JV)\{ z }) and f N+ i(zf) = «, 

and if t iiJV+ i = r* (resp. t iiJV+ i < r*) 

£(jV + 1) ={/} U B(JV) U {z} (resp. B(JV) U {2}), 

TZ(N + 1) =ft(JV) \ {z} (resp. } U K(N) \ {z}). 

Otherwise, {S*(t) : t i>N < t < r*} C »4.(iV) and S 1 * settles on a red site z' . With 
^ = /jv(^') we have < 7$ and after Sn+i(t* — £ i)iV + t^ +1 ) = S*(t*) we can set 

S N+1 (l + r* -V + ^+i) = S*,(l + U iN ), ... (3.21) 

we set t it N + i = r*, we turn blue the site z', we set f^+iiz') = i and so on. 

Since the number of red sites is finite this procedure necessarily reaches an end and one 
immediately checks that we define in this way An+i = Bn+i U7Zn+i together with the times 
{ti t N+i, i < N + 1}, and a function /jv+i : A(N + 1) — > {1, . . . , N} with all the required 
properties (/zv+i is one to one and t/ Ar+1 (*),jv+i < ^ ^ f° r an 2 m TIn+i)- 

Note first that we have the orbits inclusion by construction. Now, the law of (A(N) : 
iV > 1) is that of the internal DLA process. Indeed, the part of the flashing trajectories 
{(S*(t) : t itN <t <r*), i < N}, that can be used together with (S* N+1 (t) : < t < r£ +1 ), 
to build A(N + 1) have increments that are independent from A(N). 

Finally, for any k > 1, ipN '■ •A.ffl) — > A*(N) does associate with any site outside Uj^Sj 
a site outside Uj^Sj. Indeed, with any blue site, ipN = 9* ° In associates that site itself, 
while with each red site ip^ associates the end point of a flashing trajectory that visits that 
site. And a flashing trajectory that exits Uj<kSj necessarily settles outside Uj<fc<S,-. 

4 Estimates on the Harmonic measure 

We gather in this section two results which deal with the hitting probability of sets. The first 
one relies on a discrete mean value theorem for the Green's function. This latter theorem 
relies on Green's function estimates in [7], and Proposition IA.1I given in the Appendix. The 
second result is Proposition 13. 11 which we prove in Section H~2"l The set we wish to hit is not 
a sphere, and the proof is inspired by Lemma 5 of [8], which only gives an upper bound. 
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4.1 A discrete mean value theorem 

Our main result in this section is the following. 

Theorem 4.1 Let {A n ,n £ N} be a positive sequence with A n < Kn 1 ^ 3 for some constant 
K , and set r n = n — A n . There is a constant K a , such that for any A C dM n 

E[M(\B rn \l ,n,A))-E[M{M rn ,n,A)) | < K a \A\. (4.1) 

Written explicitly, (14. ip reads 

\\M rn \xF (S(H n )=z*)- J2 F y( S ( H n) = z*)\<K a . (4.2) 



We now recall a classical decomposition (Lemma 6.3.6 of [H]). For a finite subset A, y G A, 
and z* G dA 

By (14. 2 p and (14.31) with A = B n , we have reduced Theorem 14.11 to proving a discrete mean 
value theorem which we formulate next. We keep the same notation as Theorem 14.11 

Proposition 4.2 For z G M n , and n — ||,z|| < 1, 

| |B r J x G n (0,z) - G n(y,z) I < K a . (4.4) 

ym rn 

Remark 4.3 Note that a related (but distinct) property was also at the heart of [8]. Namely, 
for e > 0, and n large enough, if z G B n , and n — ||z|| > en, 

\M n \xG n (0,z)>J2G n (y,z). (4.5) 



Proof. We use an improved version of Lemma 2 of [7\. Using G n (0,z) = G(0,z) — 
E z [G(0, S(H n ))] (Proposition 1.5.8 of [6]), and (12.71) . one obtains by a Taylor expansion 
that for a constant K\ (independent on n) 

\ Ud G n (0,z)-2^\<^, where a(z) = E z [\\S(H n )\\ - \\z\\] . (4.6) 

Now, r d = n d — dA n n d ~ l +0(A 2 n n d ~ 2 ), so that using (14.61) . and the hypothesis A n = (^(n 1 / 3 ), 
and < n — \\z\\ < 1 



lrJG n (0,z) = (ri + O(rt 1 )) + 



nr rr 



{n d - dA n n d -' + 0(Aln d - 2 ) + O^ 1 )) ( 2^4 + °(~ d ) 

\ n n 

2a(z)(n-dA n ) + 0(1) 



(4.7) 
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A martingale argument (Lemma 3 of [7]) yields for a constant K\ 

| G n (y,z)-2a(z)n\ < K h (4.8) 



Proposition IA.1I of the Appendix reads here as follows. There is K b such that for z G B n 
with n — \\z\\ < 1 

| Yl G n (y,z)-2a (z)dA n \<K b , where a {z) =E z [\\3{H n )\\ - \\z\\\H n < H{B r J] . 

y£S(r n ,n) 

(4.9) 

Now, combining (14. 8p and (14.91) we obtain (when < n — \\z\\ < 1) 

| J2 G ™(^> z ) - 2n ( a ( z ) - a ^ z )) dA A <K t + K b . (4.10) 



Now, we combine (14.61) and (I4.10p we obtain for a constant K 2 , 

l r JG n (0,z)- J2 G n (y,z) \ +2(a (z)-a(z))dA n < K 2 . (4.11) 

V&r n J 

We now bound |ao(^) — oi{z)\ by the following expression 

P 2 (H(B rn ) < H n ) x (a (^) +E 2 [\\S(H n )\\ - \\z\\\H n > H(B rn )}) . (4.12) 
Now, it is a classical estimate (see (1A.6|) ) that there is K such that for any z G A(n — l,n), 

F z {H(B rn )<H n )<^. (4.13) 

Thus, 

A n \a (z) - a(z)\ < 2A n F z {H{B Tn ) < H n ) < 2K . (4.14) 

The desired result follows at once. I 



4.2 Proof of Proposition 13.11 

For j > 0, consider Zj in Sy. We show that for all z* in C(zj) and for suitable positive 
constants oti, a 2 , 

|<P 2j (^.) = ^)<| (4.15) 

First, z* = Zj is a flashing position when Xj = 1. This happens with probability 1 / /z^, and 
gives the result. Now, consider z* G C(zj)\{zj}. We recall that the unbiased Bernoulli Yj 
decides whether we flash on dB(zj, Rj) or on dA(rj — Rj, Tj + Rj), where Rj has density gj 
given in (13. 7p . 
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Step 1: Proof of the upper bound in (14.15ft . The following obvious facts follow from 
Lemma 12.81 

(i) z* G dM(zj, \\z* - Zj\\), and z* G dAfa - \\\z*\\ - rj\,rj + \\\z*\\ - rj\). 

(ii) z* &dB(zj, \\z* -Zj\\ - 1), and z* dA( rj - \ \\z*\\ - Tj\ + 1, r 5 + | ||z*|| - rj\ - 1). 
This means that if Yj = 1, then Rj G [\\z* — Zj\\ — 1, \\z* — Zj\\[, whereas if Yj = 0, then 
£ [I II ~ r j\ ~ 1) I Ik* II ~~ r j\[- Thus, there is a constant C such that 

(i) P(Yj = l,Rj G [\\z* - Zj \\ - 1, \\z* - Zj \\[) < C 112 *'^ - , (4.16) 

and 

(ii) =0,i?, G [| |k* || -r,\ -h\\\z*\\ -rv|[) < \ (4.17) 

In the case G dM(zj,dRj), the upper bound (14.181) then follows from (i) of (I4.16p . and 
( 12. lip of Section I2~2l We consider now Yj = 0. To simplify the notation we set for h > 0, 

V h = A(r { -h,n + h), and V h = A(r { - - , + -), 

and define two stopping times 

r = inf {n > : S(n) G 0X> h U { Zj }} , and r + = inf {n > 1 : S(n) eD^U {z,-}} . 

It is enough to prove that for some constant c, and for h such that z* G <9X\, (and /i G 
[|||^||- ri |-l,|P*||- ri |[) 

(S(H(V c h )) = z*) < j^zi- (4.18) 
We use a last exit decomposition, and the strong Markov property to get 

¥ z . (S(H(V c h )) = z*) < G Vh (z„ Zj )W z * (H(V h ) < r + ) max F x (S(r) = Zj ) 

V / x£&D c h 

= F z * (H{V h ) < r + ) max G Vh (x, Zj ) (4.19) 



< P z * (H{V h ) < t + ) max G(x, zj 



xedv c h 



It follows, from a Gambler's ruin estimate, that for a constant Kq 



P z . < r+) < ^. (4.20) 

Now, from the Greens' function asymptotics (12.71) 

sup G(x, Zj ) < sup ( ° d + j-^) ■ (4-21) 

Note that the distance between Zj and T>h is of order h. We use (14.201) and ( I4.2ip in (I4.19p 
to obtain (fl~T8l) . 

Step 2: Proof of the lower bound in (14.151) . Note the following two facts. 
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(iii) By Lemma [2.8[ z* has a nearest neighbor, say z, in M(zj, h) with 

1 



\z — Zj\\ < Ik* - Zj\\ — 



Tj\ < \\\Z*\\ - rj \ 



Ay/d 

This means that if h G [\\z* - Zj\\ - l/(4y/d), \\z* - Zj\\[, then z* G dM(zj, h). 

(iv) By Lemma 12.81 z* has a nearest neighbor, say z, in A(rj — h, rj + h) with 

1 

Wd 

This means that if h G [\\\z*\\ —rj\- l/(4\/tZ), ||2*|| - r,-|[, then z* G <92\. 

We deal separately with the cases \ \\z*\\ — Tj\ < hj/2 and |||z*|| — Tj\ > /ij/2. 
Consider first the case |||z*|| — rj\ < hj/2. On the event Yj = 1, and for h such that 
z* G 9B(zj, /i), we have 

P(y j = 1, i?, G [P*-^||-l/(4Vd), K-^||D > -73-P«, MHidBfah))) = z*) > - rT . 

Aij tlj 

(4.22) 

Thus, for some constant ot\ (that depends on d), we have (14.151) . 

Consider now the case |||z*|| —Tj\ > hj/2. It is enough to prove, for h such that z* G <92\,, 
and for some constant c (that depends on d) 

P Zj (S(H(Vl)) = z*)>-^. (4.23) 

Let be the closest site of dE(0, rj) to the segment [0,2;*], and x* be the closest site of 
a©(0,rj + /i/ 2 ) to the segment [0,z*]. We set T = B(x*, \\z* - x*\\) n XV It may be that 
M(x* , \\z* — x* ||) fl = and if so, one would only have to consider a site at a distance 1 
from z*, say z G V h , and such that B(x*, \\z — x*\\) fl = 0, and work with z instead of z* 
in the sequel. We assume henceforth that M(x*, \\z* — x*\\) fl^ = 0. 

By (14.31) with A = Z\, and the strong Markov property, 



P*. (S(H(dV h )) = z*) > G Vh (zj,Zj)F z * (H(T) < r+) minP, (S(r) 



xer 



> P 2 * (H(T) < r+) mmGv 

' x& 



h [X, Zj). 



(4.24) 



Since z* G C(zj), y* and Zj can be connected by 20 overlapping balls of radius hj/10 in 
such a way that, applying Harnack's inequality 20 times (see Theorem 6.3.9 in [§]) to the 
harmonic functions Gv h {x, •), we can estimate from below the last factor in (14.241) . There is 
a constant K2 such that 

min G Vh (x, Zj) >c§ min Gv h (x,y*) 

>c% x ™l v G(x,y')-KAG(S(H(V> h )),y>)] (4 25) 



h d ~ 2 \(^2/2) d ~ 2 ) ~ h d j 



d-2 ■ 
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As a consequence of (14.251) . we just need to prove that the first factor in (I4.24p is of order 
1/hj at least. We realize the event {H(T) < r + } in two moves: we first hit the sphere 
M(x* , R* /2), and then we exit from the cap dM(x* , R*) which lies in T>h- 

P z . (H(T) < r+) >P 2 * (H(M(x*,R*/2)) < H(M c (x*,R*))) 

«*tf* m p >( H (*< x 'W A ) = w.«*») (4 ' 26) 

We invoke again Harnack's inequality to have for y £ dB {x* , R*/2) 



F y (H (B c (x*, R*) n V h j = H (B c (x*, R*))j > c H F (H (E n B c (x*, R*)) = H (M c (x*, R*))) . 

(4.27) 

We invoke now (12.111) to obtain for some constant K 3 

P (H (B c (x*, R*) nV h )=H mx*, R*))) > ci m ^^ hl > K 3 . (4.28) 
We gather now (I4.26p . (I4.27P and (I4.28P to obtain the desired lower bound. 



5 The flashing process fluctuations 

In this section we prove Theorems 11.41 and 11.51 To do so we use the construction in terms of 
exploration waves of Section 13.21 



5.1 Tiles 

We recall that we have defined a cell of Sj in (I3.10p . as the intersection of a cone with Sj. 
We need also a smaller structure. We define, for any Zj in Sj, 

Cfa) = Sj n{xeR d : 3A > 0, 3y e B( Zj , hj/S),x = \y} . (5.1) 

As in Lemma 12 in [7], concerning locally finite coverings, we claim that, for h Q large enough, 
there exist a positive constants c±, and, for each j > 0, a subset Ej of Sj such that 



< <*JSt and S j = |J C( Zj ). (5.2) 



h- 



For any Zj G Sj, we call tile centered at Zj, the intersections of C(zj) with Sj. We denote by 
T(-2j) a iz/e centered at Zj, and by 7j the set of tiles associated with the shell Sj-. 

T, {T(.:/) : ZjGSj}. (5.3) 

Let us explain the reason for /ij/5 in the definition of a tile. It implies a fundamental feature 
of the flashing process. For any z £ <Sj, there is £ Sj such that 

*eP|{C(v): yeT(%)}. (5.4) 
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Indeed, let Zj G £j be the site realizing the minimum of {||z — y|| : y G Sj}. There is A > 
and u G B(zj, 1), such that z = Xu. Now, there is Zj G Ej such that \\zj — Zj\\ < hj/5, and 
for any y G T(%), we have — < 2hj/5. Thus 

M Gfl{^,|): yGT(%)}^,Gfl{C(y): 2/ G T(£,)} . 

5.2 The inner ball 

For n > 0, we take iV = |B n |, we recall that A*(N) = \J k >iAl(N), and write A* instead of 
A*(N). We consider 

T* = min {k > 1 : U j<k Sj A* k } . (5.5) 

We have, for Z with < n, 

P (1(0, n + ft,) £ .A*) < P ( T * =k + l) (5.6) 

k<l 

and we estimate from above the probability P(T* = k + 1) assuming < n. For > 1 
and A C Sfc, we call Wit (A) the number of unsettled explorers that stand in A after the fc-th 
wave, that is 

N 

Wfc(A) = £l A (&(i)). (5.7) 

i=l 

We now look at the crossings of tiles of Tk- On the one hand, we will use that if Wk(T) is 
large, then it is unlikely that a hole appears in the cell containing T . On the other hand, if 
r k is small it is unlikely that Wk(T) is small. We make precise what we intend by small and 
large. For this purpose, we will show in (15. 16ft that for some constant k% > 0, and any tile 

E[W(T)]>K 1 (n — r k )h k ~ 1 , and we define h = n 1 + T > sup h k . (5.8) 

k:r^<n 

For any positive constant A, we write 

P (T* = k + 1) =P (T* = k + 1, 3T G T fc , W fc (T) < Ah d k logn) 

+ P(r = Hl,VTe T fc , W*(T) > ^ log n) , 

and we estimate separately each term in the right hand side of (15.91) . 

Estimating the first term. We show here that (for K\ and h appearing in (|5.8|1 ) for k 
such that 

2A 

r k < n ft, logn, (5-9) 

there is a constant ^2 > 0, and n large enough, such that 

P(T* = k + l,BT G T k ,W k (T) < Ah d k \ogn) < \S k \ exp (-k 2 A 2 log 2 n) . (5.10) 
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On {T* = k + 1}, we have A* k = B(0,r fc - h k ). On {T* = k + 1}, and for any T C %, 
we consider a variable L k (T) = M(M(0,r k — h k ),r k ,T) independent of W k (T), and define 
M fc (r) = W k (T) + L k (T). We have the equality in law, on {T* = k + 1}, 

M k (T) l = M(Nl {0} ,r k ,T), and W fc (T) = M k (T) - L k {T). (5.11) 



As a consequence 

P(T* = Hl,3reT fcl ^(r) < log n) < |E fc | maxP (M k (T) - L k (T) < Ah d k log n) . 

(5.12) 

We first estimate the number of explorers stopped on a tile T of %. 

Now, Mj.(T) and L k (T) are dependent random variables, but both are sums of indepen- 
dent Bernoulli variables, for which Lemma |2~U1 is designed. We introduce two notations. For 
a variable X, let X = X — E[X], and let 

2x k = E [M k (T) - L k {T)) - Ah d k login). (5.13) 

Since we need x k of (15. 13)) to be positive, we will choose A and k such that 

E [M fc (T) - L k {T)} > 2Ah d h log(n), which implies x fc > log(n). (5.14) 

Then, 

P (M k (T) - L k (T) < Ah d k logn) = P (M fc (T) - Z fc (T) < -2a*) . (5.15) 
In order to estimate E[M k (T) — L k (T)}, we invoke Theorem 14.11 with n = r k , and A n = h k 

1 /3 

(the hypothesis h k = 0(r k ) holds here). We have for some positive constants n', k±, and 
for n large enough 

E[M k (T)-L k (T)} =E[M((\M n \ - |B Pk _ h J)l 0j r k ,T)] 

+ E[M(\M n ^ hk \l , r k , T)\ - E[M(M rk _ hk , r k , T)\ 
> (|B n | - |B rfc _,J) P (S(H k ) G T) - 0{h d k - 1 ) 



> 



h 



d-l 



(5.16) 



>K\{n — r k )h d 1 . (recall that n — r k > h). 



Note that the ultimate inequality in ( 15. 161) is the estimate in ( 15.81) . In view of ( 15. 161) . 
condition (15. 14)) is ensured if A and k satisfy (1 5 . 9 j) . 

Note that ( 15161) implies that E[M k {T)\ > E[L k (T)], so that ( 12TT41 of RemarkE31requires 
only an upper bound on P[Mfc(T)]. Thus, we only treat the latter quantity. 

We distinguish two cases: (i) when r k is close to n, (ii) when r k is small compared to n. 

Step 1: We assume n — h > r k > n/2. 

We set here 2x k = Ah k logn. ( 15.14)) and (15. 1 5 j) imply that 

P (M k (T) - L k (T) < Ah d k logn) < P (M k (T) < -x k ) + P (L k (T) > x k ) . (5.17) 
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To be in the CLT regime of Lemma [2.31 when dealing with the right hand side of (I5.17p . we 
need 

< x k < E[M k (T)}. (5.18) 

Let us now estimate E[M k {T)]. Note that by using (12. lip and r k > n/2, we have for positive 
constants K X) K[ 

t. \ d-l 

hi 



K [n d -± >E[M k (T)} = \M n \¥ (S(H k ) E T) 

f h \" h " (5 ' 19) 

>K x n d > K 2 n d ^ = K 2 nh d T 1 . 

Thus, E[M k (T)} > x k and we are in the Gaussian regime for (12.121) . Similarly, as in ( 15. 19ft . 
we have a lower bound 

E[M k {T)\ < K' 2 nh d k -\ (5.20) 
Thus, there is k 2 > such that for a large enough n, 

4 f A 2 h 2d loe 2 n\ 

P (M fc (T) < -x k ) < e < exp - * * = exp (-k 2 A 2 log 2 n) . (5.21) 

V ibK 2 nh k J 

As already noted, (12.141) yields 

P(L k (T) > x k ) < exp(-K 2 A 2 log 2 n). (5.22) 

Step 2: We assume r k < n/2. 

We have, using ( 15. 16ft . for n large enough 

/ i \ d—l l / h \ ^~ 1 / / h \ ^~ 1 

(5.23) 



2x k > K\n d -r%)[ — \ - Ah a k \ogn > —rf — - Ah k log n > —ri 
We define here 



K > fhi?\ d ~ 1 k' 
Xk = —n d — , (and note that x k > —nht 1 ). (5.24) 
16 \r k J 16 

As previously, we have ( 15.171) . From ( 15.191) . we have for some positive K' 2 



E[M k (T)}<K' 2 n d [^) . (5.25) 



d-l 



Now, using Lemma 12.31 for n large enough 



P (M fc (T) - L k {T) < Ahflogn) < 2exp f-^min (l, 



%K'J J (5.26) 
< exp (— k 2 A 2 log 2 n) . 

Collecting (ICTD . (IQ2]1 and (15^61) together with (I5TT2D and flOjl . we conclude that (IBTTOD 
holds. 
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Estimating the last term. The last term in the right hand side of (15.91) is bounded using 
a simple coupon-collector argument. Indeed, the event {T* = k + 1} means that there is 
one uncovered site in S k . By (15.41) . there is Zk G such that this site is a possible settling 
position of all explorers stopped in T(zk). Now, if {Wk(T(zk)) > Ahflogn}, Proposition 13. II 
tells us that the probability of not covering this site is less than (1 — a 2 / 'h d ) to the power 
Ahf log(n). In other words, 



P (T* = k + 1, VT E T k , W k (T) > Ahflogn) < \S k \ (l - ^ 



Ahf. log n 

(5.27) 

< \S k \ exp {-a 2 A\ogn} . 



Conclusion. First, choose A large enough so that 



exp (-a 2 A\ogn) < —. (5.28) 



n 2 



Recall the decomposition (15.61 and (15.91) . and assume that r/ satisfies (15. 9p . Then, ( 15. 1Q|) and 
(15.271) yield that for n large enough 

P(B(0,ri + hi ^^*(|B n |)) < |B n | (exp(-K 2 A 2 log 2 n) +exp{-a 2 Alogn)) < \. (5.29) 

The right-hand side in (15.291) is summable, and Borel-Cantelli lemma yields the inner control 
of Theorem 11.41 

5.3 The outer ball 

This section follows closely [7]. The features of the flashing process allow for some simplifica- 
tion. We keep the notation of the previous section. There, we proved that for some positive 
constant S 

P{{M(0,n-5hlogn) c A*}) = 1 -e S (n), with e\n) < +oo. 

n>l 

As consequence, the following conditional law can be seen as a slight modification of P. 

P 5 (.) = P(-\ {B(0 1 n-5h\ogn) C A*}). (5.30) 

We begin by proving that, under P s , the probability to find some k with < 2n and 
some tile T in Tk with Wk{T) larger than or equal to 2Ah d \ogn for a large enough A 
decreases faster than any power of n. First, note that, under P s , we have 

W k {T) <M k {T)-L 5 k {T), with L s k = M(M(0,n-6h]ogn),r k ,T). (5.31) 

Now, 

P S (W k {T) > 2Ah d logn) < P 5 (M fc (T) - L{{T) > 2Ah d logn) . (5.32) 



Fluctuations for internal DLA 



22 



By Theorem 14. 1[ for some positive constants K', K and for n large enough 



h 



d-l 



E [M k {T) - L{{T)] < K' (n d -(n- 6hlogn) d ) -J^ + 0(h 



d-l\ 
k ) 



< K'n d d 



Shlogn h 



d-l 



n 



d-l\ 



(5.33) 



< K'dShhfr 1 logn + OihfT 1 ) < Kh d \ogn. 
Choosing A > K, we get for n large enough so that P(B(0, n — 5h logn) C .4.*) > 1/2 

P 5 (W fc (r) > 2Ah d \ogn) < 2P (V fe (T) - E[M k (T)] > ^h d \og^j 

+2P (lI(T) - E[L s k (T)} < ~h d \ogn)j (5.34) 

As in the previous section E[M k (T)) is of order n d h d ~ l /rfr 1 , i.e, of order nh^ 1 . In addition 
E[L 5 k (T)\ is smaller than E[M k {T)\. We conclude once again by invoking Lemma 12.31 

Now, let F k denote the event that no tile T in contains more than 2Ah d log n unsettled 
explorers after the fe-th exploration wave. We denote by T k = cr(£o, • ••,&:)> an d n °te that 
F k and {1(0, n - 5h d \og(n)) C A*} are J-^-measurable. 

For any tile T G %, let z k G T, k be such that T = T{z k ) and denote by C = C{z k ). We 
are entitled, by Proposition 13.11 to use a coupon-collector estimate on the number of settled 
explorers during the k + 1-th exploration wave. On F k n {1(0, n - 5h d \og(n)) C A*}, and 
for some positive constant K\ 



E 



A* k+1 n c 



•Fk 



>\C\\l 



KJ 



W k (T) s 



> \C\ 1 - exp 



-Oi\- 



W k {T) 

K 



4 



* M Wh{T) w k{r) 



exp < —ol\- 



W k {T) 

K 



(5.35) 



> K x W k {T) inf 

x<2A log 71 



-OL\X 



X 



We now write for some positive constant K 2 



-a\x 



inf 

x<2A log n 



> 



inf 



-a\x/2A log n 



> 



2Alogn x<2A\ogn x/2A\ogn 

1 1 - e~ ai:c if 2 
inf > 



2A\ogn x<i x 



logn 



We conclude that on F k n {1(0, n - 5h d \og(n)) C A*}, 



E 



A* k+1 n C 



> #1#2 



^(T) 
logn 



(5.36) 
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Summing over all tiles we get, for a different constant K, (because of the finite, ^-independent 
overlapping between tiles), we obtain on F k H {B(0,n — 5h d log(n)) C A*} 



E 



Also, since W k (S k ) < |B(0,n)|, 



A* k+1 n S k 



> K 



W k (S k 
logn 



& [ 1 -F k n{B(0,n-Sh d \og(n))cA*}- A -k+l 1 1 ^fcj ^ A 71 ^\ P k>- 



Since P({B(0, n - 5h d log(n)) c .A*}) > 1/2, 



logn 



E s [W k (S k )} 



logn 



2n d P{F c k ). 



(5.37) 



In other words, noting that .4.£ +1 D 5^ = WfeOSfc) — W / / c+ i(5fc + i 

if 



£ d [W fc+ i(5 fc+ i)] < 1 



logn 



E d [W k {S k )] + 2n d P{F c k 



(5.38) 



By iterating (I5.38p . we obtain that for any e, E 5 [W l+elog 2 n (S l+elog 2 „)], decreases faster than 
any power of n, when / the lowest index for which r/ > n. Also, the probability (under P !) 
of seeing at least one explorer reaching the shell S l+elog 2 n is summable. Using Borel-Cantelli 
lemma, this yields the proof of Theorem 11.41 



5.4 Optimality of the fluctuation exponent 

Let time k be such that r k = n — Ah, for a large arbitrary constant A. We show that 
P(T* = k + 1) decays faster than any polynomial in n. 

On the event {T* = k + 1}, we have, after the k-th wave and for some constant K, 

|B r J = M(B rfc ,r fc ,S fc ) = M(^,r fc ,S fc ) W A (£ fe ) = |B n | - |B r J < AKn d - 1 x h. (5.39) 

This means that there exists z k G S fe such that, for some positive constant K', 

vd-l 

W k {B(z k , 3h) n S fe ) < if'n^ 1 xh — < ifV. (5.40) 

By construction, only the explorers stopped inside B(z k ,3h) can cover C(z k ) (for n large 
enough). As a consequence, we can think of a coupon-collector problem, where an album of 
size \C(z k ) | has to be filled when we collect no more than K'h d coupons. Thus, the probability 
of {T* = k + 1} is bounded from above by the probability of filling such an album, which 
is less than exp(—c(A)h d / 2 ), for some explicit constant c(A) dependent on A. This result 
is based on the following simple coupon-collector lemma (together with Proposition 13. ip . 
which we did not find in the vast literature on such problems. 



Fluctuations for internal DLA 



24 



Lemma 5.1 Consider an album of L items for which are bought independent random cou- 
pons, each of them covering one (or possibly none) of the possible L items. IfYi is the item 
associated with the i-th coupons, we assume that for positive constants ai,a2, such that for 
any j — 1, . . . , L, 

f <m = j)<f • (5-41) 
Let tl be the number of coupons needed to complete the album. Then, for any A > 0, 

P(r L < AL) < exp ( _^L± y/Z J . (5.42) 

Proof. We denote by Oi the time needed to collect the i-th distinct item after having 
collected i — 1 distinct items. The sequence cr 2 , . . . , ctl} is not independent, but if 
Qk = ■ ■ • > Yk\), and r{k) — <j\ + - — h <7fc, then for i — 1, . . . , L 

1 - ° l(L -' + 1) ) fc > P(« > k\g T(l . 1} ) > (l - ° 2(L -' + 1) ) t . (5.43) 
Indeed, calling £(i — 1) the set of the first i — 1 collected items, 

P((7; > fc|£ T (i_i)) =P ({y r (i_l) + l, . . ., Y^-ij+fc} C £(z - l)|(? r (i_i)) 

= (^(y 6 5(i-l)M)* ( 5 - 44 ) 
= (l-P(y^(i-l)|fVi)))*- 
Using (E3ID we deduce (15^31) from pi| . Now, (15^431) gives that 

oti(L-i + l) a 2 (L — i + l) 



z — ' «2 — 2 + 1 «2 

i=VL i=\J~L 

Thus, condition (15.471) holds for P > exp(2«2^4)- Finally, note that 

max ^E[a L -i\Q T {L-i-i)\, i = v 7 ^, . . . , BVZ^ < ^ 



and set 

y Ljot\ 



E[aL-i\Q T {L-i-l)\ —<?L-i 

Jii — -, r <- 1. 



(5.45) 



as well as 

Now, we look for B such that 

By/L 

E[a L _ t ] > 2AL. (5.47) 

Note that 



Fluctuations for internal DLA 



25 



For x < 1, note that e x < 1 + x + x 2 to obtain for < A < 1, by successive conditioning 



Bv 7 ! \ / 

--y/Z J \i=VL 



< e -A«avT Yl (i + ^sup^X 2 ^^]) (5-48) 

i=y/L 

<exp ^-Aa i Av / L + A 2 ^sup J E[X J 2 |^ r(L _ i _i) 



Finally, we have, using f )5.46p . 



ByfL B^fL n, 2 

2 H a L-i\yT{L-i-\ 



J2 s*pE[X?\Gr { L-i-i)] < £ ^sup ^^g^^ < 25v/Z. (5.49) 
The results follows as we optimize on A < 1 in the upper bound in (15.481) . I 



A Time spent in an annulus (By S.Blachere) 

This appendix is devoted to an asymptotic expansion of the expected time spent in an 
annulus A(r n ,n) for r n < n, when the random walk is started at some point z within the 
annulus, and before it exits the outer shell. 

Proposition A.l Consider a sequence {A n ,n £ N} with A n < Kn 1 ^ 3 for some constant 
K. Let r n = n — A n , and z £ A(r n ,n). There is a constant Kb, independent on z and n, 
such that 

| £ G n (z,y)-(2dA n a (z)-2d(n-\\z\\) 2 )\<K b ((n-\\z\\)Vl) , (A.l) 

y£A(r„,n) 

with 

a (z) = E z [\\S(H n )\\ - \\z\\\H(B c (0,n)) <H(B(0,r n ))] . 

Remark A. 2 The statement is true in dimension 2, when Green's function is replaced by 
the potential kernel 



a{x,y) = E x 



£) 1 {5(0 = x} - 1 = y} 



(A.2) 



Proof. Our strategy is to decompose a path into successive strands lying entirely in the 
annulus. The first strand is special since the starting point is any z £ A(r n ,n). The other 
strands, if any, start all on <9B(0,r n ). We estimate the time spent inside the annulus for 
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each strand. Let us remark that we make use of three facts: (i) precise asymptotics for 
Green's function, (ii) {G(0, S(n)), n G N} is a martingale, and (iii) {^(n)]) 2 — n, n E N} 
is a martingale. 

Choose z G A(r n ,n). We define the following stopping times (D i} U i} i > 0), correspond- 
ing to the i th downward and upward crossings of the sphere of radius r n . Let 9(ri) act on 
trajectories by time-translation of n-units. Let r = H(B rn ) A H n , D = U = 0, and 

Di = rl H (B rn )<H n + ool Hn<H[Brn) . 

If Di < oo, then U\ = H rn o 6(Di) + D\, whereas if D x = oo, then we set U\ = oo. We now 
proceed by induction, and assume Di, U are defined. If A = oo, then D i+ i = oo, whereas 
if Di < oo, (and necessarily Ui < oo) then 

A+i = Ui + (rl T = H {B rn ) + ool T=Hn ) o 6{Ui), and U i+1 = D i+1 + E Tn o 9{D i+1 ). 
With this notation, we can write 

oo 

£ G n (z,y) = E z [r] +^E z [ro 0(A)1a<oc] 

(A.3) 

= E z [r]+¥ z (D 1 <oo) x |(z), 

where 

OO I — 1 

l(^) = ^ A [t ° < oo] (1 - P,p j+ i = oo\Dj < oo)) . (A.4) 

i=l j=l 

Now, we compute each term of the right hand side of (1A.3I) . 

We have divided the proof in three steps. 

Step 1: First, we show that there is a positive constant K, (independent of z and n) such 
that when z G A(r n ,n), then 



|P, (A < oo) - ^ | < 1L ((n - |k||) V 1) . (A.5) 



Note that when z G B(0,n), and n — ||z|| < 1, ( 1A.5I) yields 



E z [||S(r)|| -11^111^1 = 00] , X 
A„ 1 " A? 



|P 2 (A < 00) - z[ " A J- 1 < tj, (A.6) 



Secondly, we show that for 2 G A(r n ,n), and z > 1 

| P , n \ E z [OlgTOjl - ||g(A + i)) IIId.o^XooIA < 00] 1 if 

|P 2 (A+i = oo|A < 00) ! -| < -tj. (A. 7) 

Our starting point is the classical Gambler's ruin estimate 

p m <no)= G(0,z)-E z [G(0,g(r))|A = oo] ( , 

Z[ 1 ) E z [G(0,S(r))\D 1 <oo]-E z [G(0,S(r))\D 1 = ocY { ] 
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We now expand Green's function using asymptotics (12. 7p . For this purpose, it is convenient 
to define a random variable 



X(z) = — (\\S(t)\\ 2 - \\zf) , and to set r] 



d-2 



By expressing S{t) in terms of X(z), we have 

1 1 



\\S{r)\\ d - 2 \\z\\ d - 2 



1 + 



X(z) 



Note that for any z G A(r n ,n), X(z)/\\z\\ is small. Indeed, 

X(z) (\\S(t)\\-\\z\\)(\\S(t)\\ + \\z\ 



\\z\\ \\z\\ z 
Since A n = n — r n = 0(n l l^), we have for n large enough 



(A.9) 



(A.10) 



X(z) 2(n+l)A n 8A n 

< — : — — < . and sup 



n - A r . 



n 



z£A(r„,n) \ \\ z 



\X(z) 



< 



8 3 Al 1 



n n z 



(A.ll) 



More precisely, X(z) is of order 2(||S , (r)|| — ||^||). Indeed, A^ < K'n for some K' > 0, and 
fjAIOjl yields 



X{z) = 2{\\S{t)\\-\\z\\) + 



{\\S{r)\\-\\z\\f 



\X{z)-2{\\S{t)\\-\\z\\)\<-^. 



A^ 
(A.12) 



Finally, we have a constant such that 



X(z) 



-'/ 



< 



K 



For any z ^ 0, Green's function asymptotics (12 .7p and (IA. 13[) yields 
|G(0, S(r)) - G<0, .) - <A + !L±i|Wl 

Using (IA.8P and (IA.14|) . we obtain 

E 2 [X(^)|/J 1 = oo]-C(2) + 0(i) 



< ^ 
~~ n d ' 



(A.13) 



(A.14) 



P 2 (D l < oo) 



^ [X^IA = oo] - E z [X{z)\D x < oo] +C(z) - C(z) + 0(± 



IT, (A-15) 



where 



CM = ^ 



LDi = oo 



and C(z) = ^-E z 



xHz) 



\D 1 < oo 



Using ( lA.lip . we have some rough estimates on C and C_. For any z e A(r n ,n), 



A? 



r? 



A, 



C(z) = -M = — , and 0(z) = 0^=0 — . 



A. 2 



A., 



(A.16) 



(A.17) 
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Using (1A.12j) . we have better estimates for C and C_. 

'A 



C{z) = d { ~ n -}^ ) \ 



n 



and C(z) = d {m [ y +0 



n 



(A.18) 



The rough estimates (1A.17j) allow us to derive from (1A.15|) an estimate for ¥ Z (D\ < oo), for 
any z G A{r n , n). 



F z {D 1 < oo) 



E z [\\S(r)\\-\\z\\\D 1 = oo] + o[^ 



E z [\\S(r)\\ - \\z\\\D x = oo] - E z [||5(r)|| - \\z\l\D, < oo] + O (£) 



An(l + U 



This yields f[Q|) since a (z) < 1 + (n - \\z\\) V 1 < 2(n - ||z||) V 1. 

Case where z e dB(0,r n ). 
On {Di = oo}, we have 



(A.19) 



X(z) = 2(\\S(t)\\-\\z\\) + 0(—). 



On {Z?! < oo}, we have 



X(z) = 2(\\S(r)\\-\\z\\) + 0(-). 



n 



This implies that using (1A. 18|) 



Thus, 



F z (Di = oo) 



C(z) =d^! + 0(— ), and £(*) = 0(i) 
z n ' n 



2E z [\\z\\-\\S(t)\\\D 1 <oo]+C(z) + Q(±) 
'E z [X(z)\D 1 = oo]-E g [X(z)\D 1 <oo] + C-C + 0(- 
E z [\\z\\-\\S(t)\\\D 1 <oc}+0(±) 



A n + 0(1) 
E 2 [Hzll-ll^lllD^oo] 



A., 



Ar - 



(A.20) 



(A.21) 



(A.22) 



(A.23) 



In order to obtain (1A.7jl . we write (IA.23I) on {Di < oo}, and z = S(Ui) as follows. There is 
a constant K such that on the event {Di < oo}, 



\^S(Ui) [ll) l+ i=oo] - 



A n X P 5(t/i) (Dj < OO) 



A 2 ' 



(A.24) 
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Note that (IA.23D implies that P S (y t ) (A < oo) = 1 + 0(1/ A n ), so that (IA.24j) reads as we 
integrate over {A < oo} with respect to E z 



|P Z (A+i = oo, A < oo) - 



E 2 [l A <oo (11^)11 - ||5(r)||) W^xoo] | < KF z {D t < oo 



A r 



A? 



We obtain (lA.7p as we divide both sides of HA.25|) by P Z (A < oo). 
Step 2: We show now that for any z G A(r n , n) we have 

\E Z [t] - (dA n a (z) - 2d(n - ||*||) 2 ) \ < K ((n - \\z\\) V 1) . 

When z G B n and n — ||*|| < l. (IA.26j) reads 

\E Z [r] - (dA n a (z) - 2d{n - \\z\\f) < K. 

When z G A(r n ,n), and % > 1, we show that 

E 2 [r o^) | A < oo] E 2 [(||S(A)|| - ||S(A+i)||) 



(A.25) 



(A.26) 



(A.27) 



A., 



A < ool , . . 
~\ < A2- (A-28) 



Using that {(^(n)!! 2 — n, n G N} is a martingale, and the optional sampling theorem (see 
Lemma 2 of [8]) 



AH =E 2 [||5(r)|| 2 ] 



|*|| x E 2 [X(z)} 



= \\z\\ x (E 2 [X(*)|A = oo] P 2 (A = oo) + E 2 [X(*)|A < oo] P 2 (A < oo)) . 
Thus, using flA.150 . simple algebra yields 

E 2 [r] = II^H x ((C(z) - 0(*))P 2 (A < oo) + C(zj) + 0(1). 
By recalling ( jA.18j) and (IA.5|) 



(A.29) 



(A.30) 



E 2 [r] =d (j 



|*||-r n ) 2 -(n-||*in 2 



+ 0(A n ))('^ + 0(- (n - ||:|l,V 1 



A, 



A 2 



=d(2||*|| - n - r n )« (*) + O ((n - \\z\\) V 1) 
=dA n a (z) - 2d(n - \\z\\) 2 + O ((n - ||*||) V 1) 

Note that in the case where n — ||*|| < 1, (IA.31|) yields f)A.27j) . 
Assume now that z G dB(0,r n ). From (IA.30I) . we have 

E 2 [r] = ||*|| x ((C(z) - 0(*))P 2 (A = oo) + C(z)) +0(1). 

We use (jATTjl . flA720|) and flA~2T]) to obtain 



0(1) 



(A.31) 



(A.32) 



E 2 [r] =||*|| (dA 2 +0(A n )) 



X[N|-llg(r)HlA<oo] of J_, 
A„ l A 2 / 



0(1) 



(A.33) 



=dA n E 2 [||*|| -||S(r)|||A<oo] +0(1). 
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Now, write (1A.33j) as follows. There is a constant K such that for any z G dM(0,r n ) 

(A.34) 



E z [r] E z [\\z\\ \\S(t)\\1 Di<00 ], < K 



dAl 



A n P z (A < oo) 



Al 



Note that by (TA~23]) . we have that A n P 2 (A < oo) = A n + 0(1), and - ||^(r)||l Ol<00 | < 
1, we have 



E z [t] E z [\\z\\ - \\S(t)\\1 Di<00 ], ^ K 



A r 



< 



A 2 ' 



1 dAl 

We replace z by S(Ui) in (1A.35|) under the event {A < oo} to obtain 

,E sm [r] E sm [(\\s(u t )\\ - wsfaoowmn 



(A.35) 



dAl 



W<bb1 1 < (A.36) 

n 



We multiply both sides of ( ]A.36|) by 1d 4 <ooj take the expectation on both side of (1A.36[) . and 

divide by P 2 (A < oo) to obtain flAT28]) . 

Step 3: For i > 1, we show the following bounds 

2 > 7i > 77^=, where Ti = E z [(\\S(Ui)\\ - \\S(D i+1 )\\) l A+1<0 o|A < oo] . (A.37) 

The upper bound is obvious. For the lower bound, first we restrict to {A < oo}, so 
that Ui < oo. By Lemma 12. 1[ S(Ui) has a nearest neighbor x, within B(0, r n ) such that 
\\S(U t ) || - ||s|| > l/(2y/d), and fDT37j) is immediate. 
Step 4: We show (lA.lj) using ( 1A.3j) . For p such that 1 < p < oo, let 



i-l 



a p = ^E 2 [to 0(A) I A < oo] Y[(l -P 2 (A+i = oo|A < oo)) . 



(A.38) 



i=l 



3=1 



Now, f lA.41) reads l(z) = lim^oo cr p . We establish in this step that, for some constant K, any 
integer n 

lim|l--^|<^. (A.39) 

Once we prove (1A.39|) . we have all the bounds to estimate the right hand side of ( 1A.3j) . 
Indeed, using (1A.26|I . (1A.5j) and (1A.39jl . we have 

E z [r] + P 2 (A < oo) x |(z) =dA n a (z) - 2d(n - ||^||) 2 + O ((n - ||*||) V 1) 



+ 



«o0) 
A„, 



+ 



(rc - ||z||) V 1 
A 2 " 

1 



X 



(dA 2 +0(A n )) 



= 2rfA„a (^) -2d(n- \\z\\) 2 + O ((n - \\z\\) V 1). 

(A.40) 

In order now to prove flA.391) . we introduce first some shorthand notation. For p and j 
positive integers 



ap ~ 1 dAl 



atj = P z (Dj + i = oo\Dj < oo), and (3j 



E z [to 6(Uj)\ Dj < oo] 



dAl 



(A.41) 
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With this notation (1A.7jl and (jA.28]) read as follows. 



7,- K 7- K 2K 

\aj - -^-\ < — , and \Pj - < — , so that \ aj - pj\ < -tj. 



Let us rewrite (1A.38|1 as 



p-i 

a p = a p -i - P JJ(1 - atj). 
i=i 

In order to establish flA.390 . we show by induction that 



(A.42) 



(A.43) 



Y[(l- aj )\<e p , 



with for p > 1 



3=1 



2K F 1 - T/ . , 2K 
e p = e p _i + — _[_[(!- Qj) and e x 



A 2ixv- -7/ - A 2- 

3=1 



(A.44) 



(A.45) 



Note that it is easy to estimate e p from ( 1 A. 4511 . There is a constant k,$ such that 

P fc r^T^/P / k 



e d < ^(1 + $>xp(- £ a,)) < ^ ( 1 + £ exp ( - £ 



2K 



A 



AH 



7j 

2A„ 



2^ » 
< ^/%A n 



2Kk,£ 



Now, ( 1A.44[) holds for p — 1, and assume it holds for p — 1. Then 



p— i p— l 

3=1 3=1 



Then by (IA.42[) . we have f!A.44j) with e p satisfying (1A.45|) . 
Now (I A.39[) follows as we notice that Step 3 implies that 



(A.46) 



lim 1T(1 - or,-) = 0. 



3=1 
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